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ABSTRACT 


An interpretive program for >Jhirlwind I is (iescribed 
that will accept algebraic equations, differential eoua- 
tions, etc, expressed on exowriter punched paper tape 
in ordinary mathematical notation (within certain limits 
imposed principally hy the flexowriter) as input and 
automatically provide the desired solutions* 
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lo Introduction 


Nearly everyone who has had a problem to solve on a large- 
scale digital computing machine has probably felt that it would be 
indeed convenient if one could give the machine his problem in 
ordinary mathematical language with perhaps a suggestion for a 
method of solution* We have been motivated by a synqsathy for this 
feeling in \jriting a program for l^hirlwind I that will constitute, 
we hope, a useful step in the indicated direction* The effect of 
our program is to create a computer within a computer and the pur¬ 
pose of the present report is to describe this computer and ptrovide 
a programmer's manual for its use* 

The input to the computer is in the form of a paper tape 
punched on a flexowriter, which is a specially designed typewriter 
that produces a coded pattern of holes in the tape for each charac¬ 
ter (including backspace, carriage return, etc*) that is typed. 

The problem to be solved is typed on the flexowriter in mathemati¬ 
cal form using the conventions and notations discussed below. As 
the restilting taps is fed to the computer the coded characters are 
translated by means of our program into a set of standar.l Whirl¬ 
wind instructions* After the translation is completed, the com^ 
puter carries out the indicated instructions to solve the problem, 
presenting the answers requested by the programmer in typewritten 
form* Examples of typical problems appear throughout the report. 

We have attempted to stay as close as possible to ordijaary 
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mathematical notation; however, because of the limitations of the 
flexowriter and of ottr program itself, certain typographical con¬ 
ventions most be observed. Programming rules are discussed 
throughout the report and are svimmarized in a later section. 

2. Basic Operations 

All of the lower case letters are available as variables, and 
equations of the following form are used: 

Here x stands for any lower case letter, and the blank between the 
equals sign and the comma stands for a mathematical expression in¬ 
volving lower case letters, nvimbers, plus and minus signs with 
their usual significance, parentheses, and slanting slashes / 
indicating division. Examples of such equations are: 
a = 5, 
y = -6.5a 

b = 0.0053 (a - y) /2ay, (1) 

The coaster regards the right hand side of an equation as a 
set of computational operations to be performed. The result of the 
series of computations indicated by the right hand side of a single 
equation is always a number and after the computation is finished 
this nmber is stored in the registers assigned to the letter that 
appears on the left haiui side of the equation. Equations are 
carried out in the order in which they are written. Thus, in the 
above example the value 5 is firrt stored in the space in storage 
assigned to the variable a; it is then called forth in the second 
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equation as the number stored in the registers assigned to e, 
iDultiulied by -6,3, and the result is stored in the sr,nr^, 
storage assigned to the variable y. The comouter then proceeds 
to commute the right hand side of the third eouation, making use 
of the numbers stored in the registers assigned to the symbol® 
a and y. and store the resulting number in the snace assigns ’ to 
the va,riable b« Because of this modus operandi such eouations as 


n = n + 2, 
w = - w. , 


( 2 ) 


are meaningiul. The first of these adds the number 2 to the con¬ 
tents of the storage assigned to the triable n. subtJOoii*ntly 
storing the result in the same location, and thereby cau-as tho 
ciirrent value of n to be increased by 2; the second causes 
current value of w to he replaced by iss negative. 


It should be observed that a letter veriable roust eunecr on 


the left hand side of an eouation so that it may be assigned a 
definite value before it annears on the right hand side o-e eu 
equation. This principle is illustrated in the set of equati'^''‘3 
(l) vrhere the variable a is first set to the velue y is then 
assigned a value defending or. h nuroerica.1 const'^nt. the r.urrsr* 
value of a, and so forth On the other hand a progr’’m beginning 
n = m - 1. 

(3) 

m = 1, 

vonld bg incorrect since m vroiOd have no ■nreassignel vr'ue when 
the computer attemuted to find the value of n. The meaning of 
this restriction is, of course, that methods of so''utinn must be 
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explicitly stated in general* 

Basically, the conQmter performs the elementary operations 
of addition, subtraction, multiplication and division. Certain 
conventions of a typographical or interpretational nature should 
be mentioned in connection with these operations* 

A. Parentheses 

Only standard parentheses are available typographically. No 
braces or square brackets may be used. No more than four paren¬ 
theses may be open in a single eqviation at any one time, thus, 

X = a(b + c(d + e(f + g(h +i) U + k)^)"^)), 

(4) 

is allowed while 

X = a(b + c(d + e(f + g(h + i(j + k)^))^)), 

(5) 

is not. Stated in another way, the number of left parentheses 
minus the number of right parentheses lying to the left 
of any point of the eqviation must not exceed fourj othend.se use 
of parentheses is unrestricted. 

B* Division 

A little reflection will show that an expression sUch as 
a + 2b/c + d (6) 

is mathematically ambiguous* The computer of coinrse, assigns a 
particular interpretation to such an expression. The mles by 
which it does so can be most conveniently explained as follows: 
a plus or minus sign has the effect of separating any mathemati¬ 
cal expression into two terms to the left and the right of it* 


4 



This separation with respect to a plus or minus si^ni occurrinj? 
within parentheses only applies, of course, to the total ouantity 
within the parentheses, Thus, if an expression such as that above 
is written without use of parentheses, the plus and minus sian? 
separate the expression into seuarate terms. In the case cited, 
these terms are the quantity a, the quantity 2'b divided hy c, and 
the quantity d. If in the above example, one had wished that the 
result should he the sum of a and the quantity 2 j divided hy (c+d), 
one would have had to include the quantity c+d within parentheses: 

a + 2h/(c + d) (7) 

If, on the other hand, one had wished to consider this quantity 
as the ratio of the sum a + 2h to the sum c+d, then hoth of 
these quantities should have heen included in the parentheses, 

(a + 2h)/(c + d) (8) 

If it is merely kept in mind that pluses and minuses act as sepa¬ 
rating symhols, no confusion can arise. As a further example, one 
can write an expression such as 

cd/a(h + c) + d (9) 

The computer interprets this expression as the sum of the OTiantity 
d and the ratio whose numerator is the product cd end whose de¬ 
nominator is the product of a with the quantity h + c. 

One more form of mathematical expression involving division 
is worthy of comment. The use of more than one division within 
a sin^e mathematical term is admissihle and is correctly inter¬ 
preted hy the computer. Thus a/h/c is properly interpreted as 
the product ac divided hy h. On the other hand the expression 
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a/ (b + c/d) (10) 

is interpreted as the' ratio of a to the quantity vhich is formed 
by adding to b the ratio c over d. In the latter case, the quan¬ 
tity vithin the parentheses is itself interpreted as a single 
mathematical expression to be evaluated by the computer. In 
effect, an expression contained within parentheses is isolated 
from the rest of the equation and the normal rules of multipli¬ 
cation and division apply \ri.thout regard to the remainder of the 
equation. 

C. Humerical Quantities 

Although it has been assumed implicitly in writing the above 
raatei'ial it is worth stating explicitly that the computer correct¬ 
ly intei’prets numbers in ordinary decimal notation » Ho standard 
form is assxmied for numerical quantities by the computer. Thus 
it interprets correctly such quantities as 0.012, 012,.012, 48, 

321 etc*, without ambiguity. 

D. ExTX)nent3 

i''T;er case numbers on the flexowriter appear as exponentsj 
there is an upper case minus sign but no upper case plus sign* 

The computer ^dll handle positive (unsigned) and negative inte¬ 
gral exponents. Some examples are* 
a - 5^, 

b » (a - 2)"^, 
c ■ (ab^)”^* 

d ■ a(b + c)V (o ~ a) a“^. 
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3. The ITum'ber Syatem 

Uumhera are stored and arithmetic operations are carried ont 
hy the computer in the 24, 6 numher system. In this system n'.in- 
hers appear in the machine essentially in the form a .2^ vhere a 
is a signed 24 place "binary number with l/2<a<l ora = 0 and h 
is a signed 6 place binary ntunber. Thus, it is useless to sunply 
more than 8 significant decimal figtires, although in general no 
difficulty is caused by supplying more - roundoff occurs auto¬ 
mat ically, If numbers, including, of course, numbers produced 

on 

as the result of computations, exceed Z ^ - Z^^ (or, roughly, 

13 \ 

9.2233995 * 10 ) in absolute value, an alarm will result. Any 

integer vcp to the preceding number can be written explicitly; 
e.g,, 922331 4040000000000 is a permissible number. The same 

18 

number could also be written 9.22331404.10 . It should be noted, 

however, that the computer first reads in a number as an integer, 

placing the decimal point only after the number is completely 

read in. This means that no number is allowed which; '/without 

benefit of decimal point, would be interpreted as an integer ex- 
63 39 

ceeding 2-2 . Otherwise, numbers may be written in ouite 

arbitrary decimal f onr. as has been illustrated to some extent in 
the ureceding examples, 

4 Output 

The delayer printer is used for output. The current valu~ 
of any number of letter variables be printed out by Inserting 
the instruction PRIKT followed by the desired letters, followsd 


7 



by a period* No computation may be performed during a print instruction(e*g*, PRINT 3a is not ad¬ 
missible) and comriias may, but need not be, inserted between letters* Thus, PRINT x, y* and PRINT xy* 
would both cause the current values of x and y to be printed out* The form of the print out is best 
explained ty examples* If a * 5 cvurrently the instruction PRIOT a* would cause a = +*50000000 + 01 
to be printed out* The second number; i*e*, the + 01, is interpreted as the power of 30 to be multi¬ 
plied by the first number* Thus, if •= *05 we should have a = +*50000000 - 01 and if a = - ir we 
should have a = —31415927 + 01® When a PRINT Instmiction is performed the first and last characters 
to ajjpear as output are carriage returns* If more than 6 variables are called for on a single PRINT 
instruction, 6 will be placed on each line, except possibly for the last, of course, and there will 
be no line space between lines. For example, the instruction PRINT a, b, c, d, e, f, g, h, :» vdll 
produce: 

a’s + ao b = + ** c = + *® d = + ** e = + *» f = + ** 

g=+** h=+** i=+*o 

while the instructions PRIiiT a, b, c * PRINT d, e, f* PRINT g, h, j * will produce 

a® + ®* b“ + ** c — + ** 

d® + ** e = + ** f = + ** 

g = + ** h = + ** 

where the dots stand for appropriate nmbers* 


i = + 



5* The Instructions SP« CP. SR, aiad CR 


It has been stated that equations are carried out in the 
order in vhich they are written* It is clearly desirable to be 
able to interrupt this order from time to time. We may wish to 
make the order in which the operations are performed depend on 
the result of a previous coii 5 )utation. For this purpose, or for 
piirposes of identification in general, equations may be numbered 
with dec ima l integers between 1 and 100 inclusive. Thus we may 
write 

3 X = 2, 

1 y = 0.5 (x + 2/x), 

The first equation is then identified as eqTiation number 3, the 
second as equation number 1. The equation number always precedes 
the equation to which it refers. 

The effect of the instruction SP followed by an integer from 
1 to 100 and a comma is to "transfer control" to the eqtiation 
with the specified niimber and to proceed on from that eqmtion. 

For exan^ile, SP 3, would cause equation 3 to be executed next suod 
then the equation or instruction following equation 3 in the order 
of writing. 

In the following we use the number 5 for purposes of lllus~ 
tration; it is understood that bh ^ number that has been assigned 
to an equation may be used. The effect of the instruction CP 5, 
is the same as that of SP 5, if the last computed quantity is 
negative j while CP 5 ^ is ignored if the last computed quantity is 
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positive^ 


For example, if n has the -value 3 then in the follovdng pro¬ 
gram 

e = n - 5, 

CP 5, 

X = 3 ab, 

5j = ax, 

when the instruction CP 5, is reached the equation x = 3 ab, -will 
be skipped since the quantity last computed before the CP 5, was 
e = n- 5= 3- 5= -2. Ifn had the value 5 in the above ex¬ 
ample then e = n - 5 = 0» In this case the CP may or may not 
behave like an SP since 0 nay appear as + 0 or - 0. This am¬ 
biguity need cause no difficulty in general* If, for example, it 
is desired to go thru a certain process starting with eq-uation 1 
exactly 10 times we may write 
n = 0, 

1 


n = n + 1, 
e = n - 9*5, 
CP 1, 


The quantity n is thus used as a counter, and e will then have 
the successive values - 8*5, - 7*5, - 7*5, •••» - 1*5, + *5 and we 
obtain the desired result. Further examples will be given below* 
The instruction SR 5, causes equation 5 to be execu-ted next 
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and then, after eqxiation 5» the eqnation or instruction following 
the SR 5.« instruction CR 5. "bears the same relation to SR 5 , 

that CP 5 , does to SP 5»» 

An SP, CP, SR or CE may not itself "be assigned an equation 
number. 

6. STOP 

The instruction STOP must he placed at the end of every 
program. Its effect is to stop the computer. This instruction 
may he given an equation number so that an SP or CP may he used to 
stop the cojTOutation. 

7• The Stop Character 

Plexowriter character number ^ is the "stop code" or "stop 

character". This character is used by our program inthe follow- 

1 

ing way* It may happen, for example, that a fairly large amount 
of data will be processed in several ways at several different 
times. Rather than diaplicate the long data tape onto the process¬ 
ing program tape in each instance one may simply place a stop 
character at the end of the data tape. The effect will be to stop 
the computer temporarily to allow a progreon tape to be inserted 
into the tape reader. Computation will then be resumed when the 
WWI restart button is pushed. Other sitxiations in which this 
feature may be useful will undoubtedly occur to programmers, 

8 , Ignored Characters 

Except when they occur between the digits of a number the 
following characters are ignored; spa^e, back space, tabulate 
and carriage return. Any of these characters or any combination 



of them, inserted between digits is interpreted as a multiplica¬ 
tion sign* Thus x = 23, (i*e«, no characters between the 2 and 
the 3) would set x to the value twenty-three while x = 2 3, (i*e*, 
2 followed by space followed by 3) xrould cause x to assume the 
value 6* 

The nullify character is always ignored* so that if a space 
or other erroneous character is acciden'^^^nserted between digits 

am error may be avoided by simply nullifying the character* 

/ 

It should be pointed out that superfluous characters are not 
always apparent from a print* If the combination "space, back¬ 
space" or the combination "shift to upper case, shift to lower 
case" is inserted between digits, the combination, in either case 
wotild be interpreted as a multiplication sign* Error may again 
be avoided in these cases by ntillifying the superflvious characters* 

9* The Character » Denoting 14LLt i nlication 

The lower dot of the colon has been filed off on at least one 
of the flexowriter units available at DCL* This character, or the 
colon itself, may be used to denote multiplication} e.g., 5*10^ 
has the value 50,000 as has also 5 10^ where only one or more 
spaces or tabulates is inserted between the "5" and the "10^"« 


*Except in an SP, SR, CP, OR, or STOP instruction. If a typo¬ 
graphical error occurs in these instructions, the ./hole instruc¬ 
tion must be nullified* 
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10» Examples 


Suppose it is desired to print a table of values of cosine x 
for X = 0,.1,«2,...,1 radian and that we decide to use up to 
the x^® term in the -poMex series for cosine x. This could be 
accomplished by the folloid.ng program: 

X * 0, 

1 a = 1 - x^S + xV2-3-4 - 2-3-4*5-6 

+ xV2*3*4*5*6*7*8 - x^°/2*3*4*5*6*7*8*9*10, 
PRINT X, a. 

X = X + •!, 

e = X - 1*05, 

CP 1, 

STOP 

We might prefer an iterative scheme to produce the sane 
result* The folloiiJing scheme has the advantage that the number of 
terms in the series may be adjvisted changing a single number: 

X * 0, 

1 y = 10, 

a = 1, 

2 a = 1 - a x^y (y - 1)# 
y = y - 2, 

e = 1 - y, 

CP 2, 

PRINT X, a* 

X = X + 0*1, 
a = X - 1.05, 

CP 1, 

STOP 


IS 



11, Function Subroutings 


It is clearly de8iral)le to "build into the computer subroutines 
that comoute automatically certain frequently used functions. This 
nrocess has been begun for our computer and at present the follow¬ 


ing functions are available: 

Designation 

p2 

p5 


T 

c 

P' 

P 

P 


B 


10 

,11 

,12 

13 

,14 

,15 

,16 

,17^ 

18 

,19 

20 

.21 


P 


22 


P 


23 


Function 

sqTiare root 

sine 

cosine 

tangent 

sine ^ 

cosine ^ 

tangent ^ 

secant 

cosecant 

cotangent 

absolute value 

signun 

exponential (e^) 
negative exponential (e"^) 

sinh 

cosh 

tanh 

sech 

each 

coth 

logg 

loge 
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2 

The sine of x may be obtained by writing simply F (x) and 
such an expression may be treated arithmetically in the same way 
as a variable in an equation, for example: 

X = 2(F^(y) F^(z) + ), 

would set X equal to 2('^ y sine z + j z j ). Subject to the usual 
restriction on the use of multiple parentheses (no more thrn 4 
open at a time', computation may be carried out inside the paren¬ 
theses containing the argument of a f^anction» For example, 

X = F^(x + 2/y(a + b) z”^ ), 

■vrould set x equal to cosine (x + 2 z/y(a + b) )• Further, functions 
of functions are permitted* For exan^jle, 

X = F^(y + F^(z) ), 

would set X eqiial to tan(y + cos'’ z) or 

X = F^(F^(a) ), 

woTild set X eq\ial to |a| e 

12» Additional Variables and Variable Indices 

Subscript numbers are not available on the standard flexo- 
writer so additional syiabols for variables are obtained as follows 
Any lower case letter followed by a vertical slash followed by an 
upper case integer betv/een 0 and 1023 inclusive may be used as a 
variable in precisely the same way the lower case letters alone 
have been used above. We shall call such variables “superscript 
variables.” Thus = 2n, a|^ = a|PRIKT a, xi'^, c, dl"^^ 
are legitimate expressions. Anj'’ lower case letter may be used as 
a variable index for superscript variables in the combination 
lower case letter followed by vertical slash, folloxjed by lower 



case letter* For example, when the expression xjn is encotintered, 
the current value of n is taken as the index: if n = 6 currently 
then xjn is at this time treated as x|^> 

If, in the expression x jn, n has a non-integral value greater 
than 0 . 5 , the computer rounds n off to the nearest integer. 
Negative values of n are not permitted in this connection, nor 
are values between zero and one half incltisive, with the excep¬ 
tion that the value zero is penrJ 3sible if assigned by the speci¬ 
fic equation n ® 0,* 

The roimdoff property can often be used to advantage; e.g., 
in interpolation in a set of tabular data. Suppose that we are 
given 10 values f(Xj^), fCx^),***# ^(^0^ ® function at a set 

of equally spaced points x^, X 2 >*.»>x^q. The tabular values can 
be read in by the equations 



We also require, for the purpose, the set of integers 1 to 10, 
read in as follows: 

n » 0, 

1 n = n + 1, 
k|n = n, 

e = n - 9.5, 

CP 1, 
ll 12 

setting k| = 1, k| « 2, etc. 

To find f(x), (x^-^ x<c.x^q), by linear interpolation, the 
following program is U8ed(h * the difference between Xj^ and + 1 )* 
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j = 0.5 + (x - x|^) /h 5 
m = k|j, 
n = m + 1, 
d = j +0.5-111, 
y = d fjn + (1 - d) f|m, 

giving y as the desired result. 

Printing Superscript Variables 

The general format of the print-out of superscript variables 

is the same as for ordinary variables. The instruction PRIIIT yP. 

■will caTise 



to be printed out while the instruction PRIKT b i. will cause 



to be printed if 3 is the current value of i. 

13* SP n, etc.. 

It is already possible to change conditionally the order in 
which operations are performed by using the CP instmiction. It is 
also possible to program such a change by using the following in- 
stmction: SP followed by any lower case letter followed by a 

comma* If, for exan^ile, the variable u has the value 15 when the 
instruction SP u, is reached by the control then this instruction 
will have the effect of SP 15,. The instructions CP, SPl and CR 
may similarly be used with lower case letters. Por non-integer 
values of the index, the same restrictions and round-off proper¬ 
ties hold as for variable superscript indices. 
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Superscript variables are ordinarily assigned to 1© (magnetic 
core storage) and there is room for a total of 250 such variables* 
If additional variables are desired they may be obtained in ivjo 
ways* The instruction ASSIGN followed by a lower case letter 
followed by a vertical slash followed by an upper case integer 
between 1 and 1023 inclusive followed by a comma causes space on 
the drum to be selected for variables* For exan 5 )le, ASSIGN aj^^^ 
automatically reserves space for variables a| ^ thru a inclu¬ 
sive to the drum and takes up in the storage table for variables 
in MS the space of only one ordinary superscript variable* Once 
this assignment has been made, the superscript variables in 
question may be treated as ordinary superscript variables* 

Special provision has been tiade for assigning tables of 
values to the drum* Suppose, for example, that we wish to store 
the values of a function g at the values t^, t^* ty t^ of its 

argument* Let these respective values be assigned the. symbols 

1 2 3 lA II 

g >g >g'»gi and suppose, for definiteness, that g ® 2, 

g “ 4» gp = 6, g|^ = 0* We already have one method of assign¬ 
ing these variables to the drujii 
ASSIGN g]^, 

g 1^ * 2, g|^ » 4, g] ^ “ 6, g|^ » 0, 

Exactly the same thing may be accomplished by writing 

giN = 2, 4, 6, 0, 


IS 




Any lo^jer case letter may be used in this way but only the upper 
case letter N may follow the vertical slash. Each value, includ¬ 
ing the last, must be followed by a comma. 

Once variables have been assigned to the dnm in this fashion, 
they may be treated as ordinary variables in exactly the same way 
as variables assigned to the drum (but not given initial values) 
by the instruction ASSIGH* 

The ajs type of instruction may be used in another way to 
assign variables (with their initial values) to the drum. Suppose 
that one wishes to set ai^ = 1, a|^ = 1.5, ap = 2, ap = 2.25, 
ap = 2.5, ap = 3.5, aP = 4»5,. This may be done by writing 
a|N = 1.0 ( .5) 2. (0.25) 2.5 (D 4*5, 

That is, one writes after the equals sign, the value of aP, then 
the increment in parentheses, then the value up to which one is to 
go using this increment, and so on, terminating the instruction 
with a comma. It is understood that the "a" in the above example 
may be replaced by any lower case letter but only N may follow 
the slash. 

As another example consider the little routine given in the 

discussion of interpolation in section 13 for assigning kP = 1, 

I 2 110 

k| =2,..,k| = 10. This could be replaced by 

k|N = 1 (1) 10, 

It shotdd be noted that the operations multiplication, addi¬ 
tion, etc., involving variables on the drum will in general take 
about 30 milliseconds longer per variable on the drum than the 
same operations involving IC variables. If variables on the drum 
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are used sufficiently often in a program to iiake a significant 
difference in the computer time used, and if there is space for 
the variables in MS then the variables should be assigned, or re¬ 
assigned to MS. 

For example, suppose one wishes to introduce 100 values 
fl^ thru fj^^ of a function. If they are assigned by writing 
out 



then 100 unnumbered equations, the maximum, have been used. 
Hence we might assign the values to the drum by 



and then re«^ssign them to MS (if they were to be used sufficient¬ 
ly often to make it worth the trouble) by writing something like 
n = 0, 

1 n = n + 1, 

fin = yin, 

'e = n'- 99.5, 

CP 1, 

It is difficult to say in advance just how many variables 
may be assigned to the drum for a particular problem. There is 
a maximum of 3072 (which may be increased at a later date) but 
the drxim is shared by the equations and the drum variables of a 
program so that in no case could 3072 drum variables actually be 
used. In any case, if storage is exceeded the automatic post 
mortem (described below) will inform the programmer where and how 
the trouble arose. 


20 




15» Differential Equations 

Provision has been made for the solution of ordinary differ¬ 
ential equations ty the method of Gill - a variation of the 4th 
order Runge-Kutta technique. In order to make use of this pro¬ 
vision, certain conventions must be observed* 

1) the letter t must be used for the independent variable* 

2) the letter h must be used for the increment in t* 

3) the upper case letter D is used to denote d/dt* 

Any ordinary or superscript variable (other than t and h) 
may be used for the dependent variables and for auxiliary vari¬ 
ables* A given nth order system must be reduced to a system of 
u first order differential equations* Suppose then, this reduc¬ 
tion having been done, we let y^, y 2 > •» 7^ denote the n depen¬ 
dent variables* Our system is of the form: 

dy^/dt = fj^(t, y^, y^,* • ,7^), i = l,2,*.,n. 

The way in which the solution is to be accomplished is per¬ 
haps best explained by an example* Suppose that there are two 
differential equations, that we take 2 as the initial value for t 
and that h = 0*5* Suppose further, that we wish to print out the 
solutions for t = 2, 3, 4 up to 10, that the Initial values of 
the y^ are 0, and that 

y^; y^) - 72 * ^2 ^1^ ^2^ * “ ^1» 

Our program may then be as follow: 

t = 2, 
h = 0.5, 
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1 PRIHT t, y|^, y|^. 

e = 9.7 - t, 

CP 3, 

2 D yj^ = yl^ + 1, 

D y|^ = - y|-", 

e = - e, 

CP 2, 

SP 1, 

3 STOP 

It should be observed that t is automatically increased by 
the value h upon completion of the last equation that starts vdth 
a D« Cf course, any or all of the values of t, h and the depen¬ 
dent variables may be reset during the computation* One impor¬ 
tant restriction on the form of a program should be noted. The 
carrying forward of the solution a single step of length h actu¬ 
ally involves going through the part of the program that starts 
with the first D equation and ends with the last D equation a 
total of 4 times. These 4 sub-steps take place automatically 
but it is essential that all relevant ainciDiary computation take 
place between the first and last D equations. Suppose for exam¬ 
ple that we were to change the part of the preceding program 
between equations I and 2 as follows: 
e = 9*7 t, 

CP 3, 

2 a = yP + 1, 

D y|^ = a, 
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This chan^ffe would epuepr to "be innocuous since ”a” s to 
"be appropriately reset before each entry into the differential 
eouations. However, upon coiupletion of the second H eoua^ioT: in 
each of the first 3 sub-steps, cormutation would be resumed vith 
the first D equation (and not equatior <.1 so that "a" would not 
have the current value of 3^ +1. 

Pinal''y, it has been stated that a total o^ P'^O mS super- 
scriut variables may oe used; this is not true for uroaren? in- 
volving differential eouations. If there are k eT"”^' 15a that 
begin with a D in a urogram then P‘iO - 2k MS suuerscriut variables 
may be used, 

■ 6 Post Mortems 

Certain automatic uost-mortem features have bean inoludel 
but this uart of the comuuter is essentially still in the dasim 
stage. Roughly sneaking, if a uroararo exceeds storage in any ono 
of several possible ways the computer will tell us where and how 
storage was exceeded. 

If an arithmetic alarm occurs, the programmer will be told 
the number of the equation in vrhich the alarm occurred and the 
number of the equation that ureceded the alairm. It should be 
noted here that eouations not assieped numbers by the uroevemmer 
are automatical''y assigned numbers from 101 to 200 by the com¬ 
uuter In the order in which they ere written. 

The programmer may obtain the values of any veriabTes he 
desires after an arithmetic alarm by writing the auurouriate 


PRIIJT instruction as eouation 100 



Rtiles 


1. Every equation must end with a comma. 

2* The nullify character is always ignored except when 
it occxurs between the S and P of SP(similarly for 
SR, CR, CP, STOP). 

3. Except when they occur between the digits of a num¬ 
ber or in SP etc., the following characters are 
ignored: space, back space, tabulate and carriage 
return. Any of these characters inserted between 
numbers is interpreted as a multiplication sign. 

4* Every PRIITT instruction must end \d.th a period. 

5* No more than 4 open parentheses may be witten in 
a single equation before one is closed* 

6. The instructions SP, CP, SR, CR must be terminated 
by a comma and may not be numbered * 

7» HUNT instructions may be numbered. 

8. The instruction STOP may be numbered and nay appear 
only at the end of a pi*ogram* 

9* The instruction ASSIGN must be terminated by a comma* 

10* Each assignment of a group of variables to the drum 
(by either the ASSIGN or the xlN type of instruction) 
uses up the space ordinarily occupied by one super¬ 
script v;;:riable in the table for superscript variables 
in yi.th this understanding a total of 250 super¬ 
script v£‘”iables may be assigned to MS* 

11* A maximum of 100 numbered eq^lations and 100 unnumbered 
equations may be used in any one problem* 



18» Tenroorary Rttles and Restrictions 

As of the date at vhich this report is vritten the follow¬ 
ing restrictions are in effect* It is expected that they will 
be removed in the immediate future* 

1* An SR or CR instruction may not refer to an equation 
that makes use of a ftinction subroutine* 

2* Variables may not be assigned to the drum; that is, 
neither the ASSIGN nor the xjN type of instruction 
is in operation* ' 

3* Function subroutines F thru have not been 
tested* 

All other portions of the program are believed to be oper¬ 
ating satisfactorily but exhaustive tests have not been made* 
Information relative to any troubles encountered in the use of 
this program would be appreciated by the authors* 
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